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C"| ■ Abstract. We have studied the Jaynes-Cummings-Hubbard model for a chain with 

O-J the Kerr effect (nonlinear optical effect) through fermionic approximation. We have 

observed that the Kerr effect does not cause major changes in the energy spectrum. 
However, the phase transition properties from Mott insulator to superfluid undergoes 
^ ' significant changes due to the Kerr effect. 
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1. Introduction 

In recent years, experimental success in the engineering of high-quality micro- cavities 
with interactions between photons and atoms has opened up the possibility of using 
light-matter systems as quantum simulators for ultracold atoms in optical lattices 
[U El [3]. The simplest system of radiation- matter coupling is the interaction of a two- 
level atom with a single quantized mode of an optical cavity. An effective model for such 
interactions was proposed by Jaynes and Cummings (Jaynes-Cummings model, JCM) 
50 years ago [HOE]. 

Recently, a generalization of the JCM has been proposed in order to approach 
photon-hopping between cavities in optical lattices [7]. Essentially, this generalization, 
called the Jaynes-Cummings-Hubbard Model (JCHM), describes the competition 
between the strong atom-photon coupling and the photon-hopping between cavities. In 
particular, a quantum phase diagram involving the transition from the Mott insulator 
(MI) phase to the superfluid (SF) phase has been investigated using the JCHM [H [9]. 

The JCM can be exactly solved within the framework of rotating wave 
approximation [10]. However, JCHM does not yet have an analytical solution requiring 
approximation methods. Mering et al [H] proposed an approach which treats spin 
operators as fermionic, allowing theoreticians to solve JCHM in momentum space using 
a Fourier transform over the bosonic and fermionic operators. 

Nonlinear optical effects are usually observed in an optical lattice [121 US] • Today, 
these effects have attracted significant interest due to the possibility of producing 
entangled states, an achievement particularly important to the field of quantum 
information [HI 115]. Using optical lattices, highly applicable results have been obtained 
using the Kerr effect [161 Ej- The Kerr is an atomic third order effect that occurs in 
optical cavities [18] , which is often overlooked in cavities of high quality. Generally, it 
is important to investigate the consequences of its presence. Studies involving JCM in 
the presence of the Kerr effect have been performed, including obtaining the statistics 
of photons emitted from a cavity driven by an external laser source [19j EH ETJ. In 
the present paper, we investigate the JCHM in Kerr medium to optical chain to large 
size using the theoretical approach proposed by Mering et al [UJ. We show that the 
presence of the Kerr effect in optical lattices influences the transition from MI to SF, 
increasing critical hopping. 

This paper is structured as follows. In section E we show the Jaynes- Cummings- 
Hubbard-Kerr model describing its effective Hamiltonian. We introduce the fermionic 
approximation in section |3] In section H] we expose the results. Finally, our conclusions 
are in section |5j 
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2. The Jaynes-Cummings-Hubbard-Kerr Model 

The Hamiltonian of the Jaynes-Cummings-Hubbard-Kerr Model (JCHKM) to a chain 
of L atoms is given by (ft = 1) 

H=J2 ^f a % + e 12 °t & 7 + 3^(4^" + 
i j j 

Here = a x zkia y and o^^,* are Pauli matrices, and dj (aj) is the annihilation (creation) 
operator of the light mode in the jth cavity with frequency Uj. The frequency of the 
transition energy of the atoms is denoted by e. The light-atom coupling is represented 
by g, and td is the hopping integral between dth-neighboring cavities. The first four 
terms in Hamiltonian (pQ) represent the JCHM, and the last one is the Kerr term where 
7 is the constant related to the nonlinear response of the Kerr medium [22l [23j [2U [25] . 
For a large optical chain we can write Uj = —u> = 2t((3) and td = -4, where t = 

2 

2^3, uj z and u x are the longitudinal and transversal frequencies of light, respectively, 
and u is the mean equilibrium distance between cavities [27]. 

When the hopping td — and Kerr 7 = terms vanish, the Hamiltonian ([T]) is 
decoupled into L independent JCM Hamiltonians, which has well-known eigenstates 
[26]. If 7 7^ 0, the system remains decoupled with the same eigenstates. However, when 
td 7^ 0, the cavities become coupled, increasing the complexity of the solution due to 
the fact that we cannot write the eigenstates of the whole system as a direct product 
of single-cavity eigenstates. In this situation, an appropriate approach is the fermionic 
treatment [11] followed by a mean field approximation that disregards the momentum 
transfer between photons in the Kerr term. 

3. Fermionic Approximation 

Fermionic treatment consists in replacing the spin operators with fermionic ones. This 
approach produces exact results for the JCM. However, as observed in the previous 
section, the eigenstates of JCHM are not a direct product of single-cavity JCM 
eigenstates. In this situation, the present treatment is an approximation that has 
presented results very close to other methods such as the mean field approach [TT] . 
We can write Hamiltonian ([1]) as 



H = co ^ djdj + e ^ cjc.,- + g y^(q|c,- + djc^) 



+ 12 f d 12( & } & 3+d + a] +d dj) + 7 12 afa) 
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where a + (a ) are replaced by fermionic operators (c). Assuming that the number 
of cavities is denoted by L, and then performing a Fourier transform, 



4f E e &k, Cj = -j= ]T e 27Ti ^ c fc , 



we find 



# = ^ uj k a\a k + S-X](4cfe + a*4) 

k k 

+ e H 44 + J I] 4+g4'-gMfc' ( 2 ) 
fc k,k',q 

where Wfe = 2t \^f =l c ^§=^ - C(3) . 

Due to the Kerr term, the Hamiltonian ([2]) cannot be written as a A; summation. 
Subsequently, we employ a mean field approximation to decouple this term, assuming 
that q = or that, in other words, there is no momentum exchange. Under these 
circumstances, we will see that the solution of the Hamiltonian depends on self-consistent 
equations similar to the standard mean field approach. The last term of the Hamiltonian 
(T5]) is written as 

Y H a\a\,a k a kl = jYl 4(<^4' ~ hk')a k < = 7 ( n - ^) ^ a\a k , 

kk> *-* kk' ^ ' k 

where n = ^ fc <4«fc is the photon-density operator of the system. Therefore, in the 
large chain limit (L ^> 1) we can write the Hamiltonian as 

H = [^ka\a k + ec\c k + g(a\c k + a k c\) . 

k 

in which u k = u k + 771. Because we are interested in obtaining transition points 
on the phase diagram, we employ the fact that, in the MI phase, the mean number 
of photons in the ground state, n Q , is the same for all moments, and consequently 
(n)o = j- ^2 k (a k a k )o = n o [llj - Now, by the follow replacement method, n — > n , and 
we obtain H = J2k H k , in which 

H k = u k a\a k + ec\c k + g{a\c k + a k c\), (3) 
u k =u k + 7n . (4) 

In this way, the system is decoupled and its eigenstates become direct products of the 
individual k eigenstates, allowing us to study our system through Hamiltonian (j3J). 
We assume the notation \l,m) for the state with I fermions and m photons with k 
momentum. Note that the total number of excitations, N k = a\a k + c\c k , commutes 
with H k . Therefore, we resorting to the basis {|0, n), |1, n — 1)} because the total number 
of excitations is constant, and the ground state of the Hamiltonian ([3]) is given by 

|Vo) = a\0,n) + P\l,n- 1), 

in which 



a= _ , = 
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Figure 1. Energies for typical values of e and 7 are given by ([5]). Solid lines refer to 
t/g = 0.2 while dashed ones are related to t/g = 0.001. Note that when t/g — >• 0, the 
minimum disappears. 



and 

xi = ^+(^p) 2 . 

Afterwards, it is relatively easy to find that the ground state energy of Hamiltonian 
It is easy to find that the ground state energy of the Hamiltonian ([3]) is 

El = u k n - 2 - Xk- (5) 
Now, we can compute the mean photon number in the ground state as follows 

^0 = x0ol4 S *#o) = n\a\ 2 + (n - l)\(3\ 2 = n- \(3\ 2 . (6) 

The self-consistent mean field equation is closed by inserting ()6]) in (j4]). 

When hopping, the term is non-null in the JCHM, and there is a transition between 
the MI and SF states [TT]. The equation (jSJ) allows us to investigate the influence of the 
Kerr effect in this transition. 

4. MI-SF transition 

Figure Q] shows how energy depends on k. Observe that for t — > the energy has no 
dependence on k. Consequently, the photons are uniformly distributed in momentum 
space characterizing a MI phase, which has a gap between the chemical potentials of 
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Figure 2. Mott lobes for (a) e = and (b) e = g. ashed lines refer to 7 = while 
solid ones are related to 7/5 = 0.01. 
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Figure 3. Critical hopping in terms of e to some 7 values, and n — 1 (dotted lines), 
3 (solid lines), and 6 (dashed lines). 



the particles and the hole [TT]. However, the gap is reduced as t increases until it 
becomes null at a specific critical point, t c , appearing in the transition from the MI to 
SF phases. The chemical potential of the particle /i+ = E^ 1 — Ey, while the hole's 
potential fi~ = E% — E%~ , where k' (k) is the minimum (maximum) of energy ([5]). In 
figure [H k! = L/2 and k = or k = L for any parameters n and 7. Therefore, the 
chemical potentials are calculated through the following expressions 

= E t/2 ~ E L/2i = K ~ '■ ( 7 ) 

The well-known feature of the MI-SF transition is the Mott lobe which is exhibited 
in Figure [2] for some values of n, e, and 7. The plus (minus) sign in ([7]) corresponds to 
the upper (lower) boundary of the Mott lobe. The Kerr effect produces a shift of Mott 
lobes, which is emphasized as n increases. This result is expected because the Kerr term 
in Hamiltonian ([T]) has a n 2 dependence. The lobs are also shifted through changes in 
e in agreement with the results of non-Kerr models [llj. When /z+ = /i~ the lobe is 
closed at the critical point, t c , which depends on n, 7 and e. 

The dependence of t c on e for some values of n and 7 is exhibited in figure [31 Note 
that t c decreases as e increases meaning that critical hopping is smaller as atomic-level 
spacing increases. Furthermore, the increase of n reduces t c in accordance with non-Kerr 
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model results [IT]. Additionally, t c increases with 7, and this behavior is amplified as 
large as n, which is expected due to the n 2 dependence of the Kerr term in Hamiltonian 

©• 

5. Conclusion 

We have studied properties of the MI-SF transition on a long-range-hopping JCHK 
chain through fermionic approximation followed by a mean field approach. Despite the 
fact that the Kerr effect be very small (7 <C g, see section CP and it causes negligible 
changes in the energy spectrum (see figure [1]), it significantly changes the properties 
of the transition (see figure [2] and figure [3]). A similar feature had been exhibited 
in reference (28] where it can be observed that changes in the photon hopping range 
influence slightly the energy spectrum while modifying significantly MI-SF transition 
properties. Therefore, we wish to emphasize the importance of taking into account 
properties that induce small effects in the energy spectrum when the research focus is 
MI-SF transition properties. 
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